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ABSTRACT 

X-ray variability in Active Galactic Nuclei (AGN) is commonly analysed in terms of the Power Spectral Density (PSD). The break 
observed in the power spectrum can be interpreted as a characteristic X-ray variability time scale. Here we study variability properties 
within the framework of clumpy accretion flows, in which shocks between accreting elements account for the UV and X-ray emissions. 
We derive a characteristic X-ray time scale, r x , and compare it with the measured PSD break time scale, T B . A quite good agreement 
is found in both magnitude and trend. In particular, the model dependence on black hole mass and accretion rate precisely reproduces 
the empirical relation obtained by McHardy et al. (2006). We suggest a possible physical interpretation of the break time scale and 
briefly discuss the related aspects of optical/UV variability and correlations between different wavelengths. 
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1. Introduction 

Variability in Active Galactic Nuclei (AGN) covers a wide range 
of time scales and amplitudes over the entire electromagnetic 
spectrum. The observed light curves are characterized by aperi- 
odic and featureless fluctuations, suggesting that the variability 
mechanism is of random nature. 

Early studies of X-ray variability focused on the search 
of characteristic time scales ( |McHardy| ( |2001| l and references 
therein). One of the methods used to study the temporal struc- 
ture of the variations is the Power Spectral Density (PSD) analy- 
sis. PSD results are mainly derived from high quality X-ray light 
curves provided by the Rossi X-ray Timing Explorer (RXTE) 
and XMM-Newton. The observed power spectrum is generally 
modeled by a power law of the form P v oc y", where v - I IT is 
the temporal frequency. At high frequencies (short time scales), 
the PSD presents a steep slope of a — 2, with no characteris- 
tic time scale. At low frequencies (long time scales), the PSD 
is characterized by a slope of a — 1, representing flicker noise. 
A characteristic time scale can be defined, associated with this 
bend in the power spectrum. Estimations of the PSD bend or 
break time scale are obtained by fitting broken power laws to the 
observed PSD. 

Subsequent studies showed that the characteristic break time 
scale Tg is related to the black hole mass Mbh- the larger the 
central mass, the longer the characteristic time scale, and a linear 



scaling (of the form Tb °c Mbh) has been proposed (Markowitz 
|et al.| ( |2003| ), |Papadakis| ( |2004| l). But a non negligible scatter has 
soon been noticed, with narrow-line objects not fitting into the 
general picture. Narrow-line Seyfert 1 galaxies (NLS 1) form a 
particular class of AGNs, characterized by rapid and large vari- 
ability combined with other peculiar spectral properties. NLS 1 
galaxies are usually believed to be powered by small black holes 
accreting at high rates, close to the Eddington limit. The fact 
that the break time scales are shorter in NLS 1 objects when 
compared with classical broad-line galaxies of the same mass, 
suggests that the break time scale depends on a second param- 



eter, such as accretion rate or black hole spin. (McHardy et al. 
pXH), |Uttley & McHardy] pOOB) ). Following this hypothesis, 
McHardy et al. ( 2006 ) proposed a sort of 'fundamental plane' re- 
lating variability time scale, black hole mass, and accretion rate. 
In particular, they obtained an observational scaling relationship 
in which the break time scale explicitly depends on both black 
hole mass and accretion rate. 

The PSD break time scale is expected to provide insights into 
the emission mechanisms, and several physical interpretations 
have been proposed. Various accretion disc time scales have 
been compared with the measured PSD time scales, but with no 
conclusive result (Markowit z et al.| ( 2003| l, Papadakis] p004)). 
The currently favored interpretation is based on the inner prop- 
agating fluctuation models (Lyubarskii (19971, Churazov et al. 

( ifoTnT i). 

Cooling processes responsible for the emitted radiation are 
known in considerable detail, while the process of transfer of 
gravitational energy into radiative energy is still poorly under- 
stood. In standard accretion disc models, it is provided by some 



form of viscosity parametrized by the a parameter (Shakura & 
|Syunyaev| ( |1973 1). Another likely source for heating of the ac- 
creting gas is given by dissipative processes, such as shocks. 
In previous papers (Courvoisier & Tiirler (20051, Ishibashi & 
Courvoisier (2009), hereafter Paper I), we discussed the cas- 
cades of shocks model in which radiation is emitted as a result 
of shocks between elements (clumps) forming the accretion flow. 
In this picture, clumps move with velocities determined by the 
gravitational field of the central black hole following different 
orbits depending on the different initial conditions. In the cen- 
tral regions, colliding clumps are characterized by high relative 
velocities. When such high-velocity flows converge, shocks will 
result. This process releases the bulk of the clump kinetic en- 
ergy converting it into radiation. Optically thick and optically 
thin shocks then account for the observed optical/UV and X-ray 
emissions, respectively. Characteristic time scales can be derived 
within this framework. Here we compare model time scales with 
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observed variability time scales and suggest a possible interpre- 
tation of the break seen in the power spectrum. 

The present paper is organized as follows. Some elements 
of clumpy accretion flows are recalled in Sect. 2. We derive X- 
ray variability time scales and present model results in Sect. 3; 
comparison with observational results are given in Sect. 4. We 
then discuss some aspects of optical/UV variability in Sect. 5 
and correlations between different wavelengths in Sect. 6. The 
main results are summarized in the Conclusion. 



2. Cascades of shocks in clumpy accretion flows 

In Paper I, we studied the properties of the inhomogeneous 
accretion flow formed by interacting clumps of matter, where 
shocks between elements and the subsequent evolution are at the 
origin of the radiation. 

A collision between two clumps of several solar masses, 
M c = M33 • 10 33 g, moving at the local free-fall velocity at 
a typical distance of ~100R S (parametrized in units of the 
Schwarzschild radius R = £ ■ R$ where £ = (uv ■ 100) leads to 
an optically thick shock. Following this optically thick shock the 
resulting gas cloud expands rapidly, with a fraction T] rad = 71/3- 1 
of the kinetic energy being radiated at the photosphere. The time 
it takes for the collision energy to be radiated is given by the 
expansion time 



r io 6 M 1/2 <r 1/4 s 



(i) 



We assume a spherical expansion with the photospheric ra- 
dius given by 



^„ = 3T0 15 M!M /4 cm 



(2) 



In the central regions, the expansion of the gas envelopes 
lead to interactions between envelopes originating from differ- 
ent events. Further shocks are then expected, with the expanding 
regions filling a volume ~ R\ mx . 

In Paper I, we distinguished two classes of objects, Class Q 
and Class S, according to the importance of the volume filling 
factor of the configuration relative to the typical system size, 
defined as 



R„ 
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(3) 



Class Q objects are characterized by a small filling factor 
(e <sc 1), a condition met for large black hole mass; Class S ob- 
jects are characterized by a large filling factor (e ~ 1) associated 
with a small central mass. We identified Class Q objects with 
massive and luminous quasars, while Class S objects were iden- 
tified with less luminous sources, such as Seyfert galaxies. In the 
following we focus on Class S objects, since timing analysis are 
mainly performed on local Seyfert 1 galaxies. Each class was 
further subdivided into two cases, Case A and Case B, depend- 
ing on whether the radiation time is longer (Case A) or shorter 
(Case B) than the accretion time. 

The second generation shocks occur in optically thin con- 
ditions. In these optically thin shocks, Coulomb collisions with 
ions, carrying most of the gravitational energy, energetize the 
electrons which cool through Compton radiation. The seed pho- 
tons for Compton upscattering are provided by the optical/UV 
photons emitted in the optically thick event. In a stationary sit- 
uation, the plasma temperature is set by the balance between 
heating and cooling rates. The electron temperature is then es- 
timated assuming an equilibrium between Coulomb heating and 



Compton cooling: L Cou iomb = L Co mpton- This gives a value for the 
average electron energy E e of the order of a few hundred keV, 
hence Compton cooling of the hot electrons is responsible for 
the observed X-ray emission. The resulting X-ray luminosities, 
L x , in the different classes and sub-cases have been calculated in 
Paper I. 

3. X-ray variability properties 

3.1. Characteristic time scale 

We introduce the heating time scale of the electrons given by 
T >™t = JWT~ ■ < 4 > 

V dt /heat 



At equilibrium this is equal to the electron Compton cooling 



time 



(dE\ 

\ d t Jcool 



(5) 



The above time scales are directly linked to the physical pro- 
cess responsible for the X-ray emission and can therefore be con- 
sidered as a characteristic X-ray time scale, defined as: 



T X ~ Theat ~ Tcool ■ 



(6) 



The numerical value is of the order of ~days: 

where the black hole mass is expressed in units of Mbh = M% ■ 
10 8 M Q , and the accretion rate in units of M — Mq ■ \M Q /yr. 

The characteristic X-ray time scale is proportional to the 
square of the central mass and inversely proportional to the ac- 
cretion rate. The shortest time scale is obtained in the case of a 
low mass object accreting at a high rate. We note that the depen- 
dence on mass is stronger than that on accretion rate. The de- 
pendence on black hole mass and accretion rate comes from the 
photon energy density given by the luminosity of the optically 
thick shocks, proportional to M, and the average distance of the 
optically thick shocks at the origin of the UV photons, propor- 
tional to M\ H . Thus X-ray variability seems to be determined by 
the source parameters, black hole mass and accretion rate. This 
may account for the different variability properties observed in 
different AGN classes. 



3.2. Luminosity variations 

In our picture the X-ray emission is produced by 
Comptonization of seed UV photons, emitted in the opti- 
cally thick shocks, by energetic electrons created in an optically 
thin event. The X-ray luminosity then depends on the properties 
of both optically thick and optically thin shocks through the 
photon energy density (u p h) and the electron number density 
(n e ), respectively. We therefore re-express the X-ray lumi- 
nosities leaving the explicit dependences on u p h and n e , using 
the photon energy density and the electron number density 
estimated in Paper I (eq. (13) and eq. (14), respectively). 

Here we parametrize these two quantities by appropriate val- 
ues (u p h = u y ■ 35.4erg/cm 3 and n e = rcg ■ 1.2 • 10 8 cirT 3 ), to 
obtain: 



Case A: (L x ) 2 4.9 • 10 43 E A p n Mey t? uv M\ic y ■ n % 



3„5/7„9/7 



erg/s 



(8) 



Case B: (L x ) = 2.1 - 10 43 £ 4 ^ eV ^M 8 2 «; 2/7 »^ /7 



erg/s (9) 
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We observe that the dependence of the X-ray luminosity on 
the electron number density is stronger than that on the photon 
energy density. This already suggests that the resulting X-ray lu- 
minosity and its variations are more sensitive to variations of the 
upscattering medium than those in the seed photon population. 

We can now relate variations in the X-ray luminosity to elec- 
tron density modulations. Assuming that the main dependence 
is on the electron density, we suppose a relation of the form 
L x (n e ) = K ■ nj 1 (where K is a numerical factor). We then con- 
sider small density perturbations An e around a mean value n e , 
and expanding to first order we obtain: 

AL X 9 An e 
L x 7 n e 



(10) 



Expression ( 10 1 directly relates density fluctuations to lumi- 
nosity variations. Modulations in the electron density can lead 
to significant X-ray variations of order unity, with density fluc- 
tuations being slightly amplified. Variations in the seed photon 
population would have a weaker effect. 

Density modulations typically occur on a time scale associ- 
ated with the inhomogeneous structure resulting from the suc- 
cession of shocks. We introduce the inhomogeneity time scale 
T inh ~ l/V, where / is a typical inhomogeneity size and V the 
medium velocity. This time scale is essentially determined by 
the local inhomogeneity properties and does not scale with the 
system parameters, such as the black hole mass. It can cover a 
broad range of time scales reflecting the inhomogeneities of the 
accretion flow. The slope of the luminosity variations as a func- 
tion of the temporal frequency, -j^-(v = V/l), is related to that of 
the density modulations as a function of the inhomogeneity size, 
—(I). Thus the shape of the power spectrum is given by the 
spectrum of density fluctuations. Electron density modulations 
can be transmitted, and thus observed as luminosity variations, 
only if the local inhomogeneity time is greater than the electron 
heating time (t,-,,/, > Tf, eat ). Fluctuations on time scales shorter 
than the heating time (t,-„/, < T/, efl/ ) cannot be propagated, they 
are smoothed out and effectively suppressed. This leads to a de- 
cline in variability power, which induces a break in the power 
spectrum. This is schematically illustrated in Figure 1. 

4. Break time scale: comparison with observations 

4. ) . Magnitude and functional dependence 

The Power Spectral Density (PSD) is defined as the modulus- 
squared of the Fourier transform of the light curve, in units of 
light curve variance per Hz. The typical power spectrum is char- 
acterized by a slope of a — 1, steepening to a slope of a — 2 
above a break frequency vb- A break time scale Tb = 1/vs is 
associated with the bend in the power spectrum and can be con- 
sidered as a characteristic X-ray variability time scale. 

Break time scales in AGNs have been studied by many au- 
thors and PSD measurements have been obtained for around 20 
sources (Uttley] ( f2007| l). Here we summarize a number of re- 
cent observational results. |Markowitz et al.| ( |2003[ ) analysed a 
sample of six Seyfert 1 galaxies and found a significant cor- 
relation between break time scale and black hole mass. The 
measured PSD break time scales were typically of the order 
of a few days, and the data could be fitted by a linear rela- 
tion of the form T fi (days) * M BH /10 6 - 5 M Q . The luminosity- 
time scale correlation was found to be negligible compared with 
the mass-time scale relation. A similar result was obtained in 
the case of broad-line galaxies with the break frequency de- 
creasing with increasing mass, following a relation of the form 
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Fig. 1. Sketch of a typical power spectrum in the (v • P y ) vs. 
v representation. Electron density modulations on time scales 
longer than the heating time (t,„/, > r/, efl ,) are effectively trans- 
mitted, leading to fluctuations in the emitted radiation. This cor- 
responds to the flat portion in the observed power spectrum. 
Rapid modulations on time scales shorter than the heating time 
(Tinh < T heat) are smoothed out, with the amplitude of fluctu- 
ations being considerably dampened. This corresponds to the 
steepening observed in the power spectrum. The suppression of 
fluctuations on time scales shorter than a critical value leads to a 
break in the variability power which may be associated with the 
observed PSD break. 
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10" 6 /(M/10 7 M Q )Hz ( |Papadakis|(2004| l). More re- 
cently, Marshall et al. (2008) reported a break time scale of ~34 



days in Mrk 509, in agreement with the mass-break time scale 
relationship discussed in the previous works. 

However, a significant scatter is observed in the linear black 
hole mass-break time scale correlation. Broad-line AGNs are 
consistent with a linear scaling of break time scales with cen- 
tral mass, but narrow-line galaxies are observed to lie systemat- 
ically above this relation. This implies that, for a given black 
hole mass, the break time scale should be shorter in narrow- 
line galaxies than in broad-line objects. Following this suppo- 
sition, |McHardy et al.| ( |2004"| suggested that the break time scale 
might depend on additional parameters such as accretion rate 
and/or black hole spin. The difference in break time scales ob- 
served in different AGN classes has been confirmed by |Uttley| 
& McHardyj d2005 ) who analysed the mass-time scale relation 



for all available objects with measured PSDs. They showed that 
for a given black hole mass, higher accretion rate objects indeed 
have shorter break time scales than lower accretion rate counter- 
parts. 

The main result, on which all observational works agree, is 
the existence of a scaling between characteristic time scale and 
black hole mass. The typical value of the break time scale is of 
the order of ~days for a 10 7 - 10 8 M o object. At fixed central 
mass, the break time scale is shorter in higher accretion rate sys- 
tems, such as NLS 1 galaxies. The numerical value of the X-ray 
time scale defined in eq. (|7]i thus matches the typically measured 
Tb value. This correspondence leads us to associate the charac- 
teristic X-ray time scale t x with the PSD break time scale Tb- 
The predicted scaling of t x with black hole mass is in agreement 
with the observed mass-break time scale correlation. For a given 
black hole mass t x decreases with increasing accretion rate, a 
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trend confirmed by the short time scales observed in narrow-line 
objects. 

The general relation between break time scale, black hole 
mass, and accretion rate has been quantified by |McHardy et al.| 
(2006). Following the hypothesis that the break time scale de- 
pends on both black hole mass and accretion rate, they obtained 
a scaling relationship of the form 



M 



1.12 



„0.98 



(ID 



where iHe ~ Li, i/Le, with L\, i the bolometric luminosity and Le 
the Eddington luminosity. The Eddington luminosity scales with 
black hole mass (Le °c Mbh) and assuming that the bolometric 
luminosity is proportional to the accretion rate (Lbai 00 M), the 
empirical relationship can be re-written as: 



m ¥h 



(12) 



We recall that the predicted dependence of the characteris- 
tic X-ray time scale on central mass and accretion rate is of the 
form: 



M 2 
M 



(13) 



The model dependence is therefore in excellent agreement 
with t he observational relation obtained by |McHardy et al.| 
( |20Q6l >. 

The empirical relation given in eq. (jTTJ is stated to be valid 
down to galactic black hole binary systems in which accretion is 
expected to be ruled by angular momentum dissipation in a disc. 
We will investigate whether our model results can be extended 
to this regime in a future work. 

4.2. Physical interpretations 

Several authors have proposed possible interpretations of the 
X-ray characteristic time scale by comparing different model 
time scales with the observed break time scales. In this per- 
spective, different accretion disc time scales have been compared 
with measured PSD break time scales (Markowitz et al. (2003), 
Papadakis (2004)). The orbital, thermal, and viscous time scales 
at different radial distances have been considered. Orbital time 
scales are too short; on the contrary, viscous time scales are 
much too long. The thermal time scale, at certain radii, gives the 
closest value to the measured break time scale. In such cases, 
thermal instabilities have been invoked as a possible source of 
variability. 

The characteristic variability time scale has also been dis- 
cussed in the framework of inner propagating accretion flow 
fluctuation models ( jLyubarskii] ( | 1 991) , [Churazov et al] ( [200 1 ) ) . 
In this picture, modulations produced at different outer radii in 
the accretion flow propagate inwards, until reaching the inner 
X-ray emitting region where they cause variations of the X-ray 
flux. These authors consider a geometry formed by a geomet- 
rically thin, optically thick disc surrounded by a geometrically 
thick corona extending out to large radial distances. The disc 
is assumed to be truncated at some radius from the centre. The 
characteristic X-ray time scale is then associated with the time 
scale of the modulations at the truncation radius of the disc. The 
difference in break time scale can be related to the location of 
the disc truncation radius: the smaller the truncation radius the 
shorter the characteristic time scale. 



In our model, the PSD break time scale is associated with 
the characteristic X-ray time scale, Tx- Since tx represents the 
heating/cooling time of electrons, the break time scale is directly 
linked to the physical process at the origin of the observed X-ray 
emission. Its numerical value and functional dependence are in 
agreement with the observed PSD break time scale. Fluctuations 
on time scales shorter than the heating time cannot be transmit- 
ted and the resulting suppression of rapid variations translates 
into a break in the power spectrum. The observed PSD break 
then reflects the suppression of short time scale variations. We 
are thus able to predict, at least qualitatively, the overall shape 
of the power spectrum and in particular the appearance of the 
break. 

Below the break frequency, variations are found to hold the 
same slope (P v oc v _1 , flicker noise) over several decades in fre- 
quency. For instance, the slope in NGC 405 1 remains unchanged 
for over four decades ( |McHardy et al.| p004)). The similar shape 
suggests that the same physical mechanism is responsible for the 
variations over this broad range of frequencies or equivalently 
time scales. In the inner propagating fluctuations model, this 
broad range of variability time scales is provided by the range of 
radii where modulations initially originate. Since different radial 
distances are associated with different time scales, the optically 
thin corona must be extended in the radial direction in order to 
account for the observed range of time scales. In the picture dis- 
cussed here, the broad range of time scales can be provided by 
the range of time scales associated with local inhomogeneities 
in the accretion flow. The inhomogeneity time is independent of 
the system parameters, while the heating time is determined by 
the properties of the source. In particular the latter is shorter for 
small central mass and/or high accretion rate objects. Since the 
inhomogeneity time does not scale with the source parameters, 
the condition t,„/, > Ti, ea t is easily met in such objects. Therefore 
the decline in variability (observed as a break in the power spec- 
trum) occurs on a shorter time scale, explaining the large rapid 
variability characteristic of NLS 1 galaxies. 

5. Optical/UV variability 

X-ray variability is also related to variations in the seed photon 
population. It is thus important to investigate possible sources 
of seed variations and to study the related optic al/UV variabil- 
ity. In general, optical/UV light curves are characterized by ir- 
regular and aperiodic fluctuations. This suggests that variability 
is due to random processes and stochastic approaches seem ap- 
propriate. In the framework of the so-called discrete event mod- 
els, variability is attributed to a superposition of independent and 
random events described by a Poisson distribution. The physical 
nature of the events can be represented by a variety of phenom- 
ena: supernova explosions, stellar collisions, and hot spots on the 
accretion disc surface. A characteristic variability time scale is 
given by the event duration, part of the observed UV variability 
is then directly related to the event rate. 

Characteristic variability time scales may be obtained from 
the study of time series through power spectral densities and 
structure function analysis. Collier & Peterson ( |2001 1 derived a 
characteristic optical/UV time scale of ~5-100 days in a sample 
of 10 AGNs, performing a structure function analysis. Studying 
the variability time scales in the framework of discrete-event 
models, Favre et al. ( 2005 \ showed a lack of (strong) dependence 
of the event duration on the source luminosity, hence on central 
mass. In their sample, the average luminosity of the sources cov- 
ered four orders of magnitude, while the event duration varied by 
only two orders of magnitude. Physical interpretations based on 
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accretion disc time scales, which scale linearly with black hole 
mass, can thus be excluded. In our model, the typical time scale 
for the radiation of the shock energy is given by the expansion 
time. A characteristic UV time scale can then be associated with 
this expansion time, which is of the order of T exp = \2Mj?£Jt d. 
The order of magnitude lies in the observed range and is deter- 
mined by collision parameters. It has no explicit dependence on 
the black hole mass, in agreement with observations of Favre 
et al. (2005 I. Shocks can be described by a Poisson distribution 
with each event rising and decaying on its own characteristic 
time scale given by the expansion time. These long-term varia- 
tions may be associated with large amplitude variations seen in 
the optical/UV range over hundreds of days. 

Discrete event models, in which a more luminous object is 
obtained simply by increasing the number of identical events, 
predicts a power law slope of -1/2 for the variability versus lu- 
minosity relation (cr(L) oc Zr 1 / 2 ). But this value is found to be 
incompatible with the observed trend which seems to favor a 
flatter dependence, with a slope of -0.08 (Paltani & Courvoisier 
1997)). Two possibilities are envisaged: either the variability 



process does not follow a Poisson distribution or the individual 
events are not identical, with the latter hypothesis being more 



plausible. Indeed, [Fa vre et aL] ( 2005| ) found a lack of correlation 
between event rate and object luminosity, suggesting that a larger 
total luminosity cannot only be explained by a greater number of 
identical events. This implies that the event energy cannot be the 
same for all events. 

In our picture, the average UV luminosity is given by the lu- 
minosity of a single event multiplied by the average number of 
collisions: {Luv) = Ljjy ■ (N c ). If the difference in total lumi- 
nosity is not due to the difference in the number of events, then 
single events should produce different luminosities. In Paper I, 
we have derived the UV luminosity of a single event in terms of 
the collision parameters: 



L uv z3-l0 45 m/3 Ml / 3 %ferg/ S . 



(14) 



The luminosity radiated in each shock is different, as each 
event is characterized by different collision parameters. In this 
case the difference in total luminosity is attributed to differences 
in the luminosity of individual events and not to the event rate. 
This can be reconciled with the observed flat dependence of vari- 
ability on luminosity. 

As the first shock occurs in optically thick conditions, the 
photospheric temperature can be estimated assuming blackbody 
emission. The post-shock temperature evolution has been stud- 
ied by Courv oisier & Tlirler | ( |2005] > . The predicted time delays 
between optical and UV light curves (increasing with the in- 
crease of the wavelength difference) correspond to the observed 
time lags. Moreover a very short lag between optical and UV 
is expected in our model since the optical/UV emissions share a 
common origin, arising from the same optically thick event. This 
property is in agreement with the observed quasi-simultaneity of 
optical and UV light curves (Courvoisier & Clavel ( 199 !) >)- The 
underlying common origin also seems to be supported by the 
observed equivalence of the mean optical and ultraviolet PSDs 
( |Collier & PetersonlpOOT) ). 

An additional source of variability may be related to the ge- 
ometry of the expansion following the optically thick shock. We 
have assumed a spherical expansion for simplicity, which is a 
very rough approximation. In reality, the expansion is expected 
to occur in a more irregular and asymmetric way, with different 
elements of the expanding envelopes having different speeds and 
thus providing a range of time scales. The inhomogeneous ex- 
pansion may then account for shorter time scale variations, and 



this could explain the short-term variability sometimes observed 
even in the optical/UV. 

The multiple contributions to the UV variability then affect 
the X-ray variability, since optically thin shocks arise from opti- 
cally thick ones, providing a connection between the two energy 
bands. 



6. Correlations between optical/UV and X-ray 
emissions 

The study of correlations between different energy bands are 
thought to provide information concerning emission processes 
and causal links relating different emission regions. Indeed, the 
UV and X-ray emissions are coupled through reprocessing of 
X-rays and/or upscattering of UV seed photons. Time lags with 
longer wavelengths leading shorter ones have been interpreted 
as fluctuations in the accretion disc, propagating from the outer 
optical/UV emitting regions toward the innermost X-ray emit- 
ting region, on the viscous or thermal time scale. Several multi- 
wavelength monitoring campaigns have searched for correla- 
tions and time lags between the two energy bands, in a number 
of AGNs. However, the sign of the time lag can be different from 
case to case, and even the existence of correlations is not always 
confirmed (INandra et al] ( [1998] ), |Maoz et al.] ([2002]), [Marshall 
|et al.| ( [2008 i). The global picture is still quite confusing and the 
problem seems not definitively settled yet. 

In our cascade model, the two emitting media are coupled 
through the succession of shocks, which also determines the 
sign of the time lag : optically thick shocks provide the seed 
photons that will later be upscattered by electrons heated in the 
optically thin shocks. Following the optically thick shock, UV 
photons are emitted and escape the region when the expanding 
clump reaches the size of the photosphere and becomes opti- 
cally thin, i.e. after T exp from the shock event. The expansion 
continues until neighboring envelopes overlap, leading to opti- 
cally thin shocks in which electrons gain energy on the charac- 
teristic heating time. X-rays are emitted as a result of the inter- 
action between photons, generated in the current event, and hot 
electrons, created in a previous optically thin event. Two dis- 
tinct time lags are then expected: a nearly zero lag due to the 
immediate Comptonization process, and a longer lag related to 
the temporal evolution of electrons, with the UV photons leading 
the X-rays. The first lag between UV and X-rays is very small, 
being of the order of the light travel time. In Class S, the second 
lag is expected to be of the order of the heating time, since op- 
tically thin shocks occur rapidly after the expanding sphere has 
become optically thin. A rough estimate of the order of magni- 
tude of the time lag between the two energy bands is given by 
t x ~ 5 M 2 /M d, of the order of a few days for a 1O 8 M object. 
A shorter lag is expected in sources with a smaller black hole. 

Significant correlations between X-rays and optical light 
curves with a delay consistent with zero lag have been recently 
reported in MR 225 1-178 ( |Arevalo et aT[ < |2008l » and Mrk 79 
( Breedt et al. ( 2009 1). Such correlations at very short time lags 
between X-ray and optical emissions may be interpreted as a re- 
sult of the Comptonization process of seed photons. 

Several cases with time lags of the order of ~days with 
the optical/UV leading the X-rays have also been observed 
(She mmer eTaLl pOOl) , |Arevalo et aL[ ( |2003] ), |Marshall etaE] 
(2008|). The direction of the time lag and the order of magni 



tude are thus compatible with our predicted lags. Marshall et al. 
(2008) report a time lag of 15 days in Mrk 509, with the opti- 
cal leading the X-rays. They attribute the shorter lags found in 
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other objects (of 1-2-days for less massive Seyferts) to a lower 
black hole mass. This is consistent with the above discussion, if 
the measured lags are associated with lags related to the heating 
time. 

In Class Q, where the filling factor is small, we need to take 
into account the different locations of optically thick and opti- 
cally thin shocks, with the corresponding electron travel time. 
Matter has to travel from the 100 Rs region where seed photons 
are emitted toward the innermost region where optically thin 
shocks take place. This is done on the local free-fall time, which 
is of the order of several months for a 10 9 M Q source. We would 
then expect an order of magnitude longer time lags in massive 
objects compared with Class S sources. Cross-correlation anal- 
ysis of 3C 273 show two peaks between UV and X-ray light 
curves, one at zero lag and the other with the UV leading the X- 
rays by -1.8 yr ( |Paltani et al.] ( [1998] ), |SoTdi et al.| ( p008) ). This 
suggests that the first lag corresponds to the Comptonization pro- 
cess, while the second one may be interpreted as the travel time, 
from the 100/?s region to the innermost region. Thus the overall 
picture seems to be supported by observations. However, one has 
to bear in mind that other physical processes might contribute to 
the observed properties (such as particular features possibly re- 
lated to a jet component) leading to non negligible uncertainties. 
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7. Conclusion 

We have considered characteristic time scales derived in the 
framework of clumpy accretion flows. The value of the char- 
acteristic X-ray time scale tx, associated with the electron heat- 
ing/cooling process, corresponds to the typically observed value 
of the PSD break time scale, Tb- The predicted time scale is 
shorter for small black hole mass and/or high accretion rate, cor- 
rectly reproducing the observed trend. The dependence on black 
hole mass and accretion rate we derived in eq. (f7|) remarkably 
agrees with the observational relation found by McHardy et al.| 



(2006), without additional parameters. 

In our picture, X-ray variability is attributed to both varia- 
tions in the seed photon population and density fluctuations in 
the upscattering medium, with a greater contribution from the 
latter. A break is expected in the power spectrum, since fluctua- 
tions on time scales shorter than the electron heating time are not 
observable. The associated X-ray time scale t x , directly related 
to the physical process of X-ray emission, may thus provide a 
possible interpretation of the PSD break time scale. 



References 

Arevalo, P., Papadakis, I., Kuhlbrodt, B., & Brinkmann, W. 2005, A&A, 430, 
435 

Arevalo, P., Uttley, P., Kaspi, S., et al. 2008, MNRAS, 389, 1479 
Breedt, E., Arevalo, P., McHardy, I. M., et al. 2009, MNRAS, 394, 427 
Churazov, E., Gilfanov, M., & Revnivtsev, M. 2001, MNRAS, 321, 759 
Collier, S. & Peterson, B. M. 2001, ApJ, 555, 775 
Courvoisier, T. J.-L. & Clavel, J. 1991, A&A, 248, 389 
Courvoisier, T. J.-L. & Tiirler, M. 2005, A&A, 444, 417 
Favre, P., Courvoisier, T. J.-L., & Paltani, S. 2005, A&A, 443, 451 
Ishibashi, W. & Courvoisier, T. J.-L. 2009, A&A, 495, 113 
Lyubarskii, Y. E. 1997, MNRAS, 292, 679 

Maoz, D., Markowitz, A., Edelson, R., & Nandra, K. 2002, AJ, 124, 1988 
Markowitz, A., Edelson, R., Vaughan, S., et al. 2003, ApJ, 593, 96 
Marshall, K., Ryle, W. T, & Miller, H. R. 2008, ApJ, 677, 880 
McHardy, I. M. 2001, in Astronomical Society of the Pacific Conference Series, 

Vol. 224, Probing the Physics of Active Galactic Nuclei, ed. B. M. Peterson, 

R. W. Pogge, & R. S. Polidan, 205-+ 
McHardy, I. M., Koerding, E., Knigge, C, Uttley, P., & Fender, R. P. 2006, 

Nature, 444, 730 



